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Abstract. Far infrared fine structure line data from the ISO archive have been extracted for several hundred YSOs 
and their outflows, including molecular (CO) outflows, optical jets and Herbig-Haro (HH) objects. Given the 
importance of these lines to astrophysics, their excitation and transfer ought to be investigated in detail and, at 
this stage, the reliability of the diagnostic power of the fine structure transitions of O I and C II has been examined. 
Several issues, such as the extremely small intensity ratios of the oxygen 63 fjxa to 145 /xm lines, are still awaiting 
an explanation. It is demonstrated that, in interstellar cloud conditions, the 145 pm line is prone to masing, 
but that this effect is likely an insufficient cause of the line ratio anomaly observed from cold dark clouds. Very 
optically thick emission could in principle also account for this, but would need similar, prohibitively high column 
densities and must therefore be abondoned as a viable explanation. One is left with [O i] 63 /xm self absorption by 
cold and tenuous foreground gas, as has been advocated for distant luminous sources. Recent observations with 
the submillimeter observatory Odin support this scenario also in the case of nearby dark molecular clouds. On 
the basis of this large statistical material we are led to conclude that in star forming regions, the [O I] and [C II] 
lines generally have only limited diagnostic value. 

Key words. ISM: atoms - clouds - jets and outflows - Stars: formation - Physical data and processes: Masers - 
line: formation 

1. Introduction roundi ng gas, forming; a Photo Dissociation Region 

(PDR; iTielens & Hollenbacbl llQsl . Both of these types 
In the optically thin limit, the far infrared lines of [OI] of sources &re knQwn to be brigM gourceg of infrared 

q ■ and [CII] provide valuable diagnostics for the densities atomic fine structure Une emission. It is thus instrumental 

03 ; and temperatures of the emitting regions. This is the case, tQ understand the details of the energy (and moment um) 

when the levels are excited by collisions with electrons injcction mto thc interste llar gas, if one wants to under- 

and/or atomic hydrogen in hot to warm interstellar re- gtand the emission and transfer of this radiation . In gen _ 

gions, generating major thermal energy sinks. This cool- ^ thig involveg & ^ of complcx assumptions and the 

mg radiation is so powerful that it is readily detectable need tQ gpedfy a krge number of parameterS; the precise 

over mtergalactic distances. Since these transitions are at valueg Qf which arg often nQt ^ knQwn Ag & bagic 

long wavelengths, they do not suffer from extinction in the processes may be conce aled or difficult to recover, 
dust enshrouded regions of ongoing star formation. 

The heating of these regions could reflect thc 

action of powerful shocks driv en by jets and out- 

n „ . „ , . i TT „ — , — ; — ; — , , T , — I as heat, without paying initially much attention to how 

flows irom young stellar objects (Hollcnbach & McKcc, . . , J ° , , , , 

I . .-J. . ,. , „„. , . . thc energy has been generated, delivered and dissipated. 

1989). Alternatively, _tUV photons Irom nearby young T , . „ , , ., ,. , . , . , , . , , , ,. 

. ., ii . In this nrst contribution, which is observational and statis- 

massive stars can pnotodissociate and neat sur- . , . c , . . . , „ 

tical m nature, we focus on the emission m the hire struc- 

Send offprint requests to: R. Liseau ture lines of 1 and C 11 from outflows in regions of star 

* Based on observations with ISO, an ESA project with formation. The statistics should disclose significant trends 

instruments funded by ESA Member States (especially the regarding the emission processes, comprehensible with a 

PI countries: France, Germany, the Netherlands and the United minimum of assumptions in simple models of the atomic 

Kingdom) and with the participation of ISAS and NASA. physics of these simple transitions. A more sophisticated 



For these reasons, we will analyze a clean and simple 
case, where we assume that the energy is locally available 
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analysis, possible only for a limited number of individual 
sources with known distances, will be deferred to a future 
paper. 

The paper is organized as follows: In Sect. 2, we de- 
scribe the data selection and in Sect. 3, we present the re- 
sults. These are discussed in Sect. 4, where we first address 
their reliability in contrast to mere spurious instrumental 
effects. We then examine various excitation models aimed 
at explaining the observations, initially using purely an- 
alytical methods, giving exact solutions, but further on 
also numerical approaches. Finally, in Sect. 5 we briefly 
summarize our main conclusions. 

2. Source selection and data reduction 

2.1. ISO archive 

Although the Infrared Space Observatory (ISO) was de- 
comissioned several years ago in 1998, much of the data 
collected during the mission still need proper evalua- 
tion and analysis. Thanks to the availability of and rel- 
atively easy access to the ISO data archive 1 these de- 
mands can be met in practice. Of the re trieved LWS 
data files (Long Wavelength Spectrometer; IClegg et all 
Il996|) , each contains spectral scans from roughly 45 /im 
to 200 /im for a particular position in the sky. The search 
key words (Scientific Category) were Interstellar Matter, 
Stellar Physics: bipolar outflows, circumstellar disks, star 
formation, young stellar objects, which resulted in a total 
of 552 LWS medium resolution observations in the data 
base. Spatial maps are counted as one observation, con- 
taining the observations of several positions in the sky 
(thus totalling 743). 

The reduction of these data is based on the Off Line 
Processing version 10 (OLP10), resulting in fully cali- 
brated spectral data sets. For [O i] 63 /im and 145 /im and 
[Cn]157/zm, the spectral line fluxes were extracted by 
Gaussian fitting of the spectral flux density at the line po- 
sition. The error estimates are based on half the peak-to- 
peak noise in the surrounding continuum and arc, as such, 
conservative. Saturated data were, of course, disregarded 
altogether. Limits are given as 3<r and have been obtained 
as 3 x rms x spl, where rms corresponds to the continuum 
flux density at the wavelength of the line and where the 
spectral resolution element spl = 0.29 /im for the 63 /im 
line and spl = 0.6 /im for the 145 /im and 157 /im lines, 
respectively. 

3. Results 

3.1. Comparison with other work 

A subset containing almost 100 ISO pointings was used 
to check for systematic trends in our measuring proce- 
dure. For these objects, line fluxes have been published 
and were obtained by a number of different individu- 
als, ensuring randomisation. Fig.^ shows this comparison 

1 http : //www. iso . vilspa. esa. es/ida/ index. html 
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Fig. 1. Comparison of line flux values obtained in this 
work to those found in the literature. For clarity, only 
the initial section of the distributions is shown (0 to 60 x 
10 -19 Wcm" 2 ), but the linear regression fits refer to the 
entire data sets (up to 400 x 10~ 19 Wcm" 2 ). The dotted 
line is of unit slope and meant to guide the eye only. Next 
to this line is the fit for [O i] 63 /im, followed by that for 
[O i] 145 /im and [C il] 157 /im, respectively. 



for the three fine structure lines [O i] 63 /im, [O i] 145 /im 
and [Cn]157/xm, from which it is apparent that sys- 
tematic trends are indeed present. Statistically, our line 
fluxes tend to be lower than the literature values, and this 
trend appears more pronounced for the [O i] 145 /im and 
[Cn] 157/im lines than for the [Ol]63/im line, with de- 
viations from unity by about 25% and 5%, respectively 

(Fig.[[ji. 

The reason for this can be understood as a 'maturity' 
effect, in the sense that the behaviour of the detectors 
became known better with time. The literature data re- 
fer all to OLP versions 7 or earlier, i.e. to versions in 
which the final LWS calibration had not yet been estab- 
lished. The archive data, on the other hand, are uniformly 
based on OLP 10, where the Relative Spectral Response 
Function (RSRF) is higher compared to earlier versions. 
Since the spectral flux density is inversely proportional to 
the RSRF, a lower value results for a given photo current. 
This effect is particularly pronounced in the LWS bands 
of the 145 /im and 157 /im lines. At the end of the mis- 
sion, the LWS seemed thus more sensitive. However, the 
refinement of these RSRFs remained within the absolute 
calibration accuracy of 30% anno unced already early in 
the mission ( Swinvard et all TlQOfJ) . 
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Fig. 2. Line ratio plot for the fine structure lines of O I and 
C II. Dots with error bars designate ratios based on differ- 
ential line fluxes (ON-OFF measurements), whereas the 
arrows indicate lower limits at the 3cr level. In the major- 
ity of observations, an OFF-source measurement was not 
made and in several cases, the OFFs refer to the average 
flux in the associated flows or, for non-mapped objects, 
the symbols refer merely to ON source measurements (see 
the text). 
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Fig. 3. Similar to Fig.0 but only ON-source measure- 
ments are shown, i.e. no subtraction of OFF or off data 
have been made. Limit symbols and error bars have been 
omitted to avoid crowding (see Fig.0). The horizontal 
dashed line indicates the optically thin regime for the 
[O I] 63/xm/[0 1] 145/xm line ratio above about ten. The 
two open symbols with vertical error bars refer to the av- 
erage values of data lying above and below the dividing 
line, respectively. 



3.2. Observed line ratios 

The ratios of the intensities of the fine structure transi- 
tions of neutral oxygen and ionized carbon arc plotted for 
outflows and their sources in Fig.[21 The upper bounds of 
a few x 10 of both ratios most likely reflect the limited 
dynamical range of the LWS. 

Many symbols with assigned error bars in Fig .El iden- 
tify flux differences (ON-OFF) prior to the ratioing, where 
the OFF refers to a reference position sufficiently distant 
from the source to contain only background radiation, if 
any. However, many of the data lack such a proper OFF- 
source measurement. In those cases where the sources had 
been mapped, it was generally found that the intensities 
in the map were much lower than that of the ON posi- 
tion. Consequently, an average flux of the map, excluding 
the ON position, is used instead. Finally, the majority of 
sources has only an ON measurement. In summary, the 
points in Fig.[3 represent a heterogeneous data set, with 
the potential risk that the quality could be insufficient for 
a statistical analysis. 

Selecting only one category of data will always intro- 
duce the same type of error, deminishing biasing effects 
on the sample. For example, the [O I] background is ob- 
served to be generally very weak (because of the relatively 
high excitation). Rather than leading to good data qual- 
ity, the subtraction of an OFF (or off) flux could intro- 



duce unwanted extra noise into the ON-OFF/ON-off pair. 
It might, therefore, be justified to use only ON data for 
the diagram and judge the quality a posteriori. Because 
[C n] 157 /Ltm emission is found almost everywhere in the 
sky, the corresponding [Oi] 63/im/[Cii] 157 fim will be an 
under-estimate and for this and other reasons (Sect. 4.6), 
results involving [Cn] 157 (im will only marginally be con- 
sidered in this paper. This analysis would have to rely on 
self-consistent model calculations and will thus be deferred 
to a future paper. 

In Fig.Ol line ratios based on only ON data 
are shown. Obviously, the two distributions shown in 
Fig. 12 and |21 respectively, are statistically indistinguish- 
able. For values of the ratio [O i] 63 /im/[0 1] 145 fim > 
10, the emissio n is likely optically thin (see be- 
low and, e.g., iTielens fc Hollenbachl Il985lt . Taking 
[Oi] 63/xm/[Oi] 145/im = 10 as a dividing line yields the 
averages (2.7 ±3.0, 15.2 ±4.4) above and (1.1 ±1.3, 5.3 ± 
2.6) below, respectively. These averages are shown as open 
symbols with vertical bars in Fig.[2]and they will find their 
application in Sect. 4.5. 

As expected (and demonstrated), the ratio 
[O i] 63/^m/[0 1] 145 (im remains essentially unaffected, 
when one neglects the OFF source subtractions. 
Therefore, the major observational result of this paper, 
namely that merely 35% of the data have ratios > 10, 
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rests on a solid foundation and forms the basis for the 
following discussion. 



4. Discussion 

4.1. Instrumental effects 
4.1.1. Detector responsivities 

Of interest here is of course the intercalibration among 
different detectors. All three lines fall onto GeGa photo- 
conductors, with the [Oi]145/im and [Cn]157/im lines 
falling onto stressed detectors. The relative trends of the 
RSRFs likely reflect also the behaviour of these detec- 
tors in an absolute sense and any instrumental origin for 
small oxygen line ratios could be safely disregarded: rela- 
tive to [O i] 63 /im, the [O i] 145 (im line would get assigned 
a lower flux than it should have, which would lead to an 
over-estimate of the Inn line ratio, in contrast to the 
observed behaviour. 

As can be seen in Fig-EI the [O i] 63 /im/[0 1] 145/im 
and [O i] 63/im/[C n] 157 /im data have different lower 
cut-offs, with the latter filling the region between 0.1 
and 1, whereas the former do not. The [Oi] 145/im and 
[Cn] 157 (im lines both fall onto the same detector with 
similar RSRFs. The cause of the differences of Fig.|2l is 
thus not instrumental in nature. We will return to this 
point in Sect. 4.6. 



4.1.2. LWS beam sizes 

Estimated sizes of the non-diffractio n-limited LWS be am 
appear non-uniform with wavelength IGtv et allEool . It 
is conceivable, therefore, that the observed line ratios in 
Figs.[5]and|21could be spuriously affected by this. However, 
in both cases, the effect would be less than a factor of 1.5, 
which would not be sufficient to explain any line ratio 
anomaly in excess of one order of magnitude. 

4.2. Theoretical line ratios 

Theoretical estimates of line intensity ratios can be ob- 
tained from excitation calculations, i.e. from the so- 
lutions to the statistical equilibrium equations (see 
the Appendix). Comparing these to the observed 
[O i] 63 /im/[0 1] 145 /im line ratios indicates that optically 
thin emission seems essentially ruled out for most temper- 
atures and densities (see the dotted curves in Fig.0}. 

The observations could be reconciled with optically 
thick emission, as indicated by the dashed curve in the 
figure. Temperatures would have to be low though, below 
100 K. There are at least two difficulties with this scenario. 

Firstly, the locus for the optically thick ratios has 
been computed for both lines being optically very thick 
(the line source functions approach the Planck functions). 
Not unexpectedly, at the low temperatures characteristic 
of dark clouds, the 145/im line would require excessive 
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Fig. 4. [O i] 63 /im/[0 1] 145 /im line intensity ratios for col- 
lisions with neutral hydrogen, H I, as the function of the 
temperature are shown by the dotted lines. The logarithms 
of the densities (in cm -3 ) are indicated above each curve. 
Optically thin ratios are given for the volume densities 



77(H) 



10 3 , 10 4 and lO^m" 3 . The broken line out- 



lines the ratio of optically thick lines of equal velocity 
width, each of which is approximated by the integrated 
Planck function at the given temperature, i.e.J^I^ dv = 
IS.vByiT). The optically thin intensity ratios for the in- 
verted levels, due to collisions with H2, are shown by the 
solid lines, with the logarithmic densities in cm -3 below. 
Observed line ratios fall roughly into the area inside the 
thin horizontal lines (compare with Fig.[3J). 



column densities, hardly ever encountered in molecular 
clouds (Fig.^J). 

Secondly, detailed and self-consistent models of ei- 
ther shocks or PDRs, having 'acceptable' total column 
densities, do not predict such low values of the FIR 
[O I] l ine ratios over the en tire feasible parameter space 
(e.g., iKaufman et"aTlll999lh An extensive discussion can 
be found in lLiseau et alJ 1)1999^ . 



4.3. Inversion of the level populations 

The line ratios in Fig.0] have been co mputed for oxygen 
atom s colliding with hydrogen atoms l|Launav fc Rouefll . 
Il977|) . At the low temperatures of the optically thick lo- 
cus, hydrogen is expected to be mostly in molecular form, 
however. In this case, collisional excitation of the [O I] lev- 
els will be with H2 (and He). The 3 Pj ground state of O I 
has energies Ej, such that E2 < E\ < Eg, each with 
statistical weight gj = 2 J + 1. The 63 /im line emission 
originates from transitions E\ — > E2, while the 145 /im 
line is due to Eq — ► E\. 
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Fig. 5. Hydrogen column densities iV(H) (in cm -2 ) ver- 
sus volume densities 71(H) (in cm -3 ) for unit line center 
optical depth in the [Oi] 145 /im line, To = 1. Collisional 
excitation with atomic hydrogen is assumed, in addition 
to a cosmic oxygen abundance and a Gaussian line width 
of lkms -1 . At temperatures below 100 K, column densi- 
ties for optically thick [O i] 145 /im emission become pro- 
hibitively large. 

Co llisional transitions J = — 1 are forbidden to first 
order ijMonteiro fc Floweilll"987lh which together with the 
fact that the Einstein transition probability Aqi < A\2 
leads to the possibilit y of inverted level populations, nj 
(see also lSmithLll969li . In this case, the inversion param- 
eter is 



go 

9i 



no 
m 



< 



(1) 



i.e. would become negative at low densities (n co ii <C 
n cr ; t ) and turn over to positive values toward the high 
density thermalization regime (« co n 3> n CT it), where the 
critical density is defined as n cr it = An/goi^) with the 
denominator being the collision rate coefficient. 

In Fig-El the inversion parameter for the oxygen levels 
J = and 1 is shown as a function of the density for colli- 
sions with HI and H2, respectively, and for three values of 
the temperature. Whereas 30/51 — no/ni > always for 
collisions with H I, this parameter is seen to be negative 
for H2 collisional excitation at temperatures exceeding a 
few times ten degrees and at densities lower than the crit- 
ical density. Optically thin line intensity ratios arc then 
expected to become different too. 

Th e collisional rates were adopted from ljaauet et all 
and computed for a thermally weighted distribu- 
tion of ortho- and para-H2, implying preferential collisions 
with para-H2 at low temperatures. For T < 10 3 K, colli- 
sions with ortho-H2 are about two times faster. Hence, in 
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Fig. 6. The level inversion parameter go/gi — no/ni for 
O I is shown for collisional excitation with atomic hydro- 
gen (dotted lines) and molecular hydrogen (solid lines), 
respectively. For collisions with H2 , this parameter is neg- 
ative at temperatures above 50 K and densities below 
10 5 cm -3 , implying inversion of the upper level popula- 
tion. At low temperatures, this parameter becomes inde- 
pendent of the density and the two curves become indis- 
tiguishable, approaching the ratio of the statistical weights 
(50/31 — no/ n i — ► 90/91 = 1/3, see the curves for 10 K). 



an ortho-dominatcd gas, the magnitude of the level inver- 
sion would become somewhat reduced. 

Helium makes up about 20% of the molecular 
gas by number, but including collisions with He 



^ " 1 to „„„ - 

(Mo nteiro fc Flowerl Il987|i alters the OI results only 
marginally. Similar is true for the He contribution to the 
excitation in an H I dominated gas. 

4.3.1. Comparison with the observations 

Fig. 01 shows [Oi] 63/xm/[Oi] 145 /im as a function of the 
temperature for collisional excitation with either HI or 
H2 (the dotted and solid curves, respectively). As ex- 
pected, collisions with H2 yield smaller line ratios, in closer 
agreement with the observations. At densities n(H2) < 
10 5 cm -3 , the loci of these optically thin ratios are cross- 
ing that of the optically thick case. 

This dependence on the exact character of the colli- 
sional excitation process reduces the diagnostic value of 
the oxygen lines. For example, a line ratio [Oi]63/im/ 
[O 1] 145 /im ~ 20 could indicate either molecular gas at 
100 K or atomic gas at 2 000K (although the line inten- 
sities would be very different). Another degeneracy exists 
for a ratio of about ten, with either optically thin or thick 
emission as apparently viable, but contradictory, options. 
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4.3.2. Masing atomic transitions 

Level inversion implies potentially laser or maser activity. 
The opacity at the line center Ao is given by 

«(vo) = '"I'll A A ° A 01 a(0) n(H) h (f - ^) (2) 

where a Gaussian profile with full width at half max- 
imum Avf w hm and center velocity vo has been assumed; 
a(O) is the abundance of oxygen relative to hydrogen and 
the other symbols have their usual meaning. With the 
inversion parameter being negative, the optical depth in 
the line becomes negative, r = k(vq) Z < 0, where Z is 
the physical length of the column along the line of sight. 
For amplified background radiation, / = I e ~ r °, a mas- 
ing 145 /im line could lead to even lower line ratios than 
indicated in Fig. [3] and potentially explain the very low 
observed values. 

The smallest possible column density iV(H) = n(H) Z 
for a gain by a factor of e 1 , i.e. for To = —1 at line center, 
is obtained for a minimal, i.e. thermal, line width and this 
case is shown in Fig.[7| for a few typical cloud densities. 
At low temperatures, characteristic of dark clouds, these 
columns are clearly in excess of what could be found in 
such objects. There, masing seems thus not capable of 
explaining any hypothetical 145 /im excess. 

This situation may change in warmer giant molecular 
clouds (FigJZJ, where masing in principle could contribute 
to the line intensity at a significant rate. In fact, Tielens 
& Hollenbach (unpublished) did notice this possibility in 
their models of the Orion PDR. We shall return to this 
point in a future paper. 

4.4. Optically thick emission 

To explore the behaviour of the lines in the optically thick 
regime we need to resort to numerical calculations. To 
keep with the statistical character of this paper, we sim- 
plify things and exploit the Sobolev approximation us- 
ing a Large Velocity Gradient (LVG) code. The basic pa- 
rameter is the column density per radial velocity interval, 
N(0)/Sv. Once specified, the emergent intensity as the 
function of the gas kinetic temperature and density is cal- 
culated (diffuse background radiation due to, e.g., dust 
can also be included). 

Results for the average density of 3 x 10 4 cm -3 are 
shown in Fig.JS] On the 100" scale of the ISO observations, 
this density is already on the high side. As the tempera- 
ture increases, all curves of Fig.|S] tend toward the locus 
of the optically very thick ratio, which should be com- 
pared to Fig. 01 (r — > oo). The recovery of this analytical 
result assures that the code behaves reliably also at these 
extremely high optical depths. 

From the figure it is evident, that in order to reduce the 
[O i] 63 [im/[0 1] 145 /im line ratio to below 10, column den- 
sities exceeding N(Ol) = 10 21 cm~ 2 would be required. 
Expressed as H column density or dust extinction, this 




100 1000 
Temperature, T (K) 



Fig. 7. The hydrogen column density, AT(H), of a parcel 
of molecular cloud gas along the line of sight to attain an 
optical depth of To = —1 in a thermally broadened 145 /im 
maser line. Logarithmic densities in cm -3 are given next 
to each curve. 

would mean AT(H) > 10 24 cm~ 2 or Ay S> 10 2 mag, rarely 
ever encountered on 100" scales. 

Fig-IHl also provides an explanation for the noticed 
paucity of data points for [Oi] 63/im/[Oi] 145 /jm< 1 in 
Fig-El Highly excessive column densities would be needed 
to populate this part of the diagram. However, taken to- 
gether, the evidence speaks not in favour of the hypothesis 
of optically thick emission to explain the observed small 
oxygen line ratios. We consider next the case of 63 /im self 
absorption due to a cold, relatively tenuous cloud compo- 
nent between the source and the observer. 

4.5. Foreground absorption 

At interstellar cloud temperatures, essentially all oxygen 
atoms reside in the lowest level of their ground state, with 
J = 2. In terms of the center opacity to of the J = 1 — 2 
[O i] 63 /im line the column density of all atoms can then 
be approximated by 

N(0 1) = 2 x 10 17 t Au fwhm cm" 2 (3) 

where the line width is expressed in kms -1 . Parameter 
values of unity and an oxygen abundance a(O) < a(O)© = 
8.5 x 10~ 4 , then correspond to a hydrogen column of 
roughly iV(H) > 2 x 10 20 cm -2 or, expressed differently, a 
visual extinction of Ay < 1 mag for oxygen depletion fac- 
tors less than five. 

These numbers indicate that reasonably small amounts 
of cold foreground gas could potentially absorb in the 
63 /im line, leaving the 145 /im line essentially unaffected 
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Fig. 8. Numerical results for oxygen line ratios at char- 
acteristic interstellar cloud temperatures. Displayed are 
the column densities of atomic oxygen N(0), for 5v = 
1 kms , needed to produce the line ratios for an H2 den- 
sity of 3 x 10 4 cm" 3 . For these parameters and at 20 K, line 
center optical depths in this plot range from minro = 0.05 
to max To = 5 x 10 5 in the 145 /xm and in the 63 (im lines, 
respectively. At 100 K, these values are 26 and in excess 
of 4 x 10 5 . 



of water, H 2 (Tm — lm l. lLarsson et alJ l)2005|) found that 
an equivalent column density of iV(H2) = 6 x 10 20 cm" 2 
of 10 K gas would be sufficient to account for the zero flux 
of the narrow self absorption in the H ?Q line. For an oxy- 
gen a bundance of a(O) = 3.0 x 10~ 4 ijSavaee fc SembachL 
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cm 



ll996|) . this translates to an iV(OI) = 1 x 10 
Comparing with the numerical constant in Eq. 3, this leads 
to an estimate of the optical depth of ro = 2 for the width 
of 1.1 kms" 1 . This would be sufficient to explain also the 
observed oxygen line ratios and we take this as indepen- 
dent support for the idea of foreground absorption of the 
63 fjja line. 

However, the fact that the analysis of 
[0 1] 63/^m/[0 1] 145 (im data may require detailed 
modeling implies that their usefulness as diagnostics 
of the physical conditions in the emitting regions is 
significantly reduced. 



4.6. The [Oi] 63 fim/[Cii] 157 fim intensity ratio 

We have shown that, for many sources, it is not the 145 fim 
line which is relatively too strong, but most likely the 
63 jttm line which is too weak. 

This would of course also affect the 
[Oi] 63/mi/[Cii] 157/mi line ratio and could in part 
be the reason for the apparent distortion of the distribu- 
tion to lower values (see Fig.EJ). 

An immediate implication would be that also the 
[Oi] 63/mi/[Cii] 157/xm ratio would have its diagnostic 
power largely reduced. 



and thereby altering the line ratios. Using the average val- 
ues of the [Oi] 63^m/[Oi] 145 fim ratios (Sect. 3.2) would 
indicate that a reduction of the 63 fim line by a factor of 
three (2.9 ± 1.6) would be required to reconcile the ob- 
served low ratios with the average optically thin emission. 
This corresponds to average optical depths of order one to 
two (To = l.Oio'g). The profile of any cold foreground gas 
cloud is likely to be narrow, corresponding to Auf w hm ~ 
1 to 2 kms -1 , say, whereas the line of the warmer back- 
ground is likely broader (more turbulence, high velocity 
outflow gas), providin g a quasi-continuum fo r the absorb- 
ing material (see, e.g. lPoglitsch et allll996j) . Apart from 
being an ad hoc component, this seems to work fine. 

Conceptual difficulties with this scenario could arise 
for nearby dark clouds. These are cold, exhibiting 
barely any temperature contrast, and the lines are 
generally narrow. This should conceivably result in 
very little absorption, yet, for objects of this type 
small [O i] 63 fim/\Q 1] 145 fim ratios h ave been observed 
l|Caux et all Il99a iLiseau et allll999j) . 

However, from detailed modeling of the observation 
with Odin 2 of a nearby dark cloud in the ground state line 

2 Odin is a Swedish-led satellite project funded jointly by the 
Swedish National Space Board (SNSB), the Canadian Space 
Agency (CSA), the National Technology Agency of Finland 
(Tekes) and Centre National d'Etude Spatiale (CNES). 



5. Conclusions 

Based on the examination of an extensive observational 
ISO-LWS material of [O I] and [C II] emission from out- 
flows in star forming regions, we conclude the following: 

• A large majority of observed intensity ratios 
[0 1] 63 fim/[0 1] 145 /im appears incompatible with op- 
tically thin emission, but could be reconciled with op- 
tically very thick line emission at relatively low tem- 
peratures. 

• Collisional excitation of the inverted levels of the 
ground state of O I by H2 can lead to inversion of the 
level populations. We show, however, that a masing 
145 /im line alone is likely not sufficient to explain the 
line ratio anomaly observed in dark clouds. 

• Likewise, optically thick emssion in the lines would re- 
quire unfeasibly large column densities to be consistent 
with the observations. 

• We provide evidence that foreground absorption, as 
has previously been advocated for the observed [O I] 
emission from distant luminous objects, is a viable 
explanation for observed oxygen line ratios also from 
dark clouds. 

• In galactic regions of outflows from young stellar ob- 
jects, the fine structure lines of O I and C II have gen- 
erally only limited diagnostic value. 
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Appendix A: Statistical equilibrium 

In the steady state (dmjdt — 0), and using common notation, 
the statistical equilibrium equations can be written in terms of 
the fractional populat i on nu mbers fj — nj/n to t, as (see, e.g., 
also IWatson fc Stor"el 



- /] fk "-con qkj+y^ f>< = 

k k>j 



Poglitsch A., Herrmann F., Genzel R., et al., 1996, ApJ, 462, 
L43 

Savage B.D. & Sembach K.R., 1996, ARA&A, 34, 279 
Schoier F.L. van der Tak F.F.S., van Dishoeck E.F. & Black 

J.H., 2005, A&A, 432 
Smith H.A., 1969, ApJ, 158, 371 

Swinyard B.M., Clegg P.E., Ade P.A.R., et al., 1996, A&A, 
315, L43 

Tielens A.G.G.M. & Hollenbach D., 1985, ApJ, 291, 722 
Watson D.M. & Storey J.W.V., 1980, International Journal of 
Infrared and Millimeter Waves, Vol. 1(4), 609 



k k<j 



For the 3 Pj inverted ground state of O I with energies E2 < 
Ei < Et), the exact analytical solution of this three level system 
is 



fo = 



(a + be) d 
ec-g 



fl = {a+b i e -h 

ec — b 

/2 = l-(/ +/l) 

where 
a = n co n Q2i 

b — ?i co ii qio (n co n goi — A01) 
c = 7i co ii (qio + qi2) + A12 
d = n co ii qio 

e = n co n (q i + 902) + A i + A 02 

g = n C oi\ qio (n co ii 901 + -4oi) 
q-2t) 



h = 



qio 



(A.2) 

(A.3) 
(A.4) 

(A.5) 
(A.6) 
(A.7) 
(A.8) 
(A.9) 
(A.10) 
(A.11) 



Recently, t he relevant atomic data have been 
compiled by ISchoier et alj (I2005T) in a convenient 
form and made readily accessible via the internet at 
http : //www . strw . leidenuniv . nl/ ^moldata/ 
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